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Abstract 

We study the fluctuations and correlations of conserved charges, such as the baryon number, 
the electric charge and the strangeness, at the finite temperature and nonzero baryon chemical 
potential in an effective model. The fluctuations are calculated up to the fourth-order and the 
correlations to the third-order. We find that the second-order fluctuations and correlations have 
a peak or valley structure when the chiral phase transition takes place with the increase of the 
baryon chemical potential; the third-order fluctuations and correlations change their signs during 
the chiral phase transition and the fourth-order fluctuations have two maximum and one minimum. 
We also depict contour plots of various fluctuations and correlations of conserved charges in the 
plane of temperature and baryon chemical potential. It is found that higher order fluctuations and 
correlations of conserved charges are superior to the second-order ones to be used to search for the 
critical point in heavy ion collision experiments. 
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I. INTRODUCTION 



Studies of QCD thermodynamics and QCD phase diagram have attracted lots of atten- 

ormed 



tions in recent years. People believe that deconfined quark gluon plasma (QGP) is 
in ultrarelativistic heavy ion collisions |l]-8j]. Various field theory models studies j9-|l7j in 



dicate that there is a critical point in the QCD phase diagram in the plane of temperature 
and baryon chemical potential (l8j . which separates the first-order phase transition at high 
baryon chemical potential from the continuous crossover at high temperature. Although 
there is no definite evidence that the QCD critical point also exists in the lattice QCD 
calculations, due to the sign problem at finite chemical potential, some lattice groups find 
that the QCD critical point maybe exist in the phase diagram, based on extrapolating re- 
sults at small fi/T (ratio of the chemical potential and the temperature) to those at large 
fj,/T I9M2II]. In the meantime, experiments with the goal to search for the QCD criti- 
cal point are planed and underway at the Relativistic Heavy Ion Collider (RHIC) at the 



Brookhaven Nationa 
CERN in Geneva 22 



Laboratory (BNL) and at the Super Proton Synchrotron (SPS) at 



24] 



In order to map the QCD phase diagram, locating the QCD critical point is a crucial 
and vital task. It has been known that the event-by-event fluctuations of various particle 
multiplicities are enhanced in heavy ion collisions that freeze out near the critical point. 
Therefore the QCD critical point can be found through the non-monotonic behavior of 
various fluctuation observables as a function of varying control parameters 25M34I] . Partic- 



ularly, fluctuations of conserved charges, such as the baryon number, electric charge, and 
strangeness, deserve more attentions. On the one hand, the fluctuations of conserved charges 
are sensitive to the structure of the thermal strongly interac ting matter and behave differ- 



ently between the hadronic and QGP phases 
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30 



35J. On the other hand, since 



the conserved charges are conserved through the evolution of the fire ball, the fluctuations 
of conserved charges can be measured in heavy ion collision experiments. 

In our previous work 36|, we have studied the fluctuations and correlations of conserved 
charges in the 2+1 flavor Polyakov-Nambu-Jona-Lasinio (PNJL) model at finite tempera- 
ture and zero chemical potential. We made an interesting comparison with the recent lattice 
calculations which were performed with an improved staggered fermion action at two values 
of the lattice cutoff with almost physical up and down quark masses and a physical value for 
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the strange quark mass [37J . It has been seen that our calculated results are well consistent 
with those obtained in lattice calculations, which indicates that the 2+1 flavor PNJL model 
is well applicable to study the fluctuations and correlations of conserved charges. The valid- 
ity of this effective model is expected, since the critical behavior of the QCD phase transition 
is governed by the universality class of the chiral symmetry, which is kept in this model. 
Furthermore, compared with the conventional Nambu-Jona-Lasinio model, the PNJL model 
not only has the chiral symmetry and its dynamical breaking mechanism, but also includes 



the effect of color confinement through the Polyakov loop 
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46|. 



In this work, we will extend our previous work 36] to the cases with nonzero baryon 
chemical potential. Such an extension is quite nontrivial, which allows us to use the fluc- 
tuations and correlations of conserved charges to explore the critical behavior of the QCD 
critical point and to find out the location of the critical point. Furthermore, when the 
baryon chemical potential does not vanish, the fluctuations and correlations of odd order 
develop finite values, which are identical to zero at vanishing baryon chemical potential. In 
this paper we shall pay our most attentions on studies of the non-monotonic behavior of the 
fluctuations and correlations of conserved charges near the QCD critical point and depict 
the contour plots of various fluctuations and correlations in the plane of temperature and 
baryon chemical potential. 

This paper is organized as follows. In Sec. II we introduce the fluctuations and correla- 
tions of conserved charges and simply review the formalism of the 2+1 flavor PNJL model. 
In Sec. Ill we give our calculated numerical results of the fluctuations of conserved charges. 
In Sec. IV we present the numerical results of the correlations among conserved charges. 
Our summary and conclusions are given in Sec. V. 

II. FLUCTUATIONS AND CORRELATIONS 

In this section we focus on the cumulants of the conserved charge multiplicity distribu- 
tions, which can be expressed as the derivative of the pressure (P) of a thermodynamical 
system with respective to its chemical potentials corresponding to conserved charges, i.e. 

BQS = d^+\P/T±) 

Xi * h ~ d^ B /Tyd^ Q /Tyd^ s /Tr 1 ] 
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where T is the temperature, and Hb,q,s are the chemical potentials for baryon number, 
electric charge, and strangeness, respectively. They are related with the quark chemical 
potentials through the following relations, 

12 11 ,11 

[in = -Hb + g//Q, Hd = - -Hq, and /x fl = -/m b - —fj,Q - ft s . (2) 

where fi u ,d,s are the chemical potentials for u, d, and s quarks, respectively. Denoting the 
ensemble average of conserved charge number Nx {X = B, Q, S) with (Nx), we can obtain 
the second and higher order fluctuations of the conserved charges as follow 

= ^Nx 2 ), (3) 
xf = ^(^x 3 ), (4) 
xf = ^((^x 4 )-3(5iV x 2 ) 2 ), (5) 

where 5Nx = Nx — (Nx) and V is the volume of the system. In the same way, we can 
also obtain mixed cumulants of the conserved charge distributions, which are known as the 
correlations among conserved charges. For example, 

Xfi y = ^pWxSNy), (6) 

X?/ = ^- 3 (SNx6N Y "), (7) 
X^V = ^(5N X 5N Y 5N Z ), (8) 

It should be noted that when the chemical potentials are vanishing, i.e. Hb,q,s = 0> the 
fluctuations and correlations of conserved charges in Eq.([T]) (also known as generalized sus- 
ceptibilities) are nonvanishing only when i + j + k is even, but when /!# has finite values, 
which is the case discussed in this paper, the generalized susceptibilities also develop finite 
values when i + j + k is odd. 

We now adopt the 2+1 flavor Polyakov-loop improved NJL model to study the fluctua- 
tions and correlations of conserved charges near the QCD critical point, which is a nontrivial 



extension to our previous work 36(, where the fluctuations and correlations of conserved 



charges were calculated in the 2+1 flavor PNJL model at finite temperature but with van 



ishing chemical potentials. It is interesting to notice that the calculated results in 36[ are 



well consistent with those obtained in lattice calculations 



371 ]. which shows that the 2+1 



flavor PNJL model is well applicable to study the cumulants of conserved charge multiplicity 



distributions. Before a detailed study, let us give a brief review on the 2+1 flavor PNJL 



model for completeness. Details about the model can be found in Ref. 17 1. 
The Lagrangian density for the 2+1 flavor PNJL model is given as 

8 

£ PN jl = i ) ( i ln D ' 1 + 7oA - rho)4> + G ^2 [(^"aVO 2 + ($ij5T a il>y 

a=0 

-K [det/^l + 75)^) + det / (^(l - 75)^)] " $* , T), 
where ip = (ip u , ipd, ips) T is the three-flavor quark field, and 
D^ = d^-iA" with A» = 6{!A , A 



-iA. 



4- 



(9) 



(10) 



where A a 's are the Gell-Mann matrices in color space and g is the gauge coupling strength. 
m = diag(mQ, m,Q, rag) is the three-flavor current quark mass matrix. Throughout this 
work, we take rag = rriQ = m , while keep m,Q being larger than ra , which breaks the 
SU(3)f symmetry. The matrix fi = diag(// u , fid, fa) denotes the quark chemical potentials 
which are related with the conserved charge chemical potentials through relations in Eq.(j2J). 

In the above PNJL Lagrangian, T) is the Polyakov-loop effective potential, 

which is expressed in terms of the traced Polyakov-loop $ = (Tr c L)/N c and its conjugate 
$* = (Tr c L*) / N c with the Polyakov-loop L being a matrix in color space given explicitly by 



L(x) = Pexp 



dr AJx, t) 



l 1 u i s±4 




exp[2/3A 4 ] 



(11) 



with (3 = 1/T being the inverse of temperature and Aa_ = iA . 

In our work, we use the Polyakov-loop effective potential which is a polynomial in $ and 



4l|, given by 



W($,$*,T) 6 2 (T) 



6 4 



rp4 



with 



6 2 (T) = a + ai(^) + a 2 (^) + a 3 (^) . 



(12) 



(13) 



The parameters in the effective potential are fitted to reproduce the thermodynamical be- 
havior of the pure-gauge QCD obtained from the lattice simulations, and their values are 
given in Table [H The parameter T is the critical temperature for the deconfinement phase 
transition to take place in pure-gauge QCD and To is chosen to be 270 MeV according to 
the lattice calculations. 
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TABLE I: Parameters for the Polyakov-loop effective potential U 



Of) 




n 2 


0,3 


63 


64 


6.75 


-1.95 


2.625 


-7.44 


0.75 


7.5 



In the mean field approximation, the thermodynamical potential density for the 2+1 
flavor quark system is given by 




+ _ i n [1 + 3$* e -( E P-^/)/T + 3$ e -2(£/-/^)/T + e -3(4- W )/Tj 
3 

+ Z l n [1 + 3$ e -(-^+A*/)/T + 3$* e -2(B|+^)/T + e -3(S/+ M/ )/Tj I 

3 J 

+2G(0 U 2 + 0/ + S 2 ) - 4JT0 U 0, S + $*, T), (14) 

where 0j's (2 = w, d, s) are the quark chiral condensates, and the energy-momentum 
dispersion relation is E l p = \J p 2 + M 2 , with the constituent mass being 

Mi = ml - 4G& + 2K<f )j <j> k . (15) 

Minimizing the thermodynamical potential in Eq. ( ll~4"|) with respective to <f) u , <frd, <j>si 
and $*, we obtain a set of equations for the minimal conditions, which can be solved as 
functions of temperature T and three conserved charge chemical potentials /zg, Hq, and /is- 

We will use the method of Taylor expansion to compute the fluctuations and correlations 
of conserved charges in the PNJL model. Before numerical calculations, we should fix the 
five parameters in the quark sector of the model, whose values used usually in the literatures 
are those obtained in Ref. fl, m l = 5.5 MeV, m s = 140.7 MeV, GA 2 = 1.835, KA 5 = 
12.36, and A = 602.3 MeV, which are fixed by fitting the observables m n = 135.0 MeV, 
m K = 497.7 MeV, Try = 957.8 MeV, and f n = 92.4 MeV. 

III. NUMERICAL RESULTS OF FLUCTUATIONS OF CONSERVED CHARGES 

In this section, we are going to present the numerical results for the conserved charge 
fluctuations in the 2+1 flavor PNJL model. It is found that the QCD critical point is lo- 
cated at about T c = 160 MeV and /xb c = 819 MeV (/zq = fis — 0) with input parameters 
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FIG. 1: (color online). Quadratic (top), cubic (middle), and quartic (bottom) fluctuations of 
baryon number as functions of the baryon chemical potential [Ib (mq = us = 0) with several 
values of temperature in the PNJL model. 

given above. This QCD critical point separates the first-order chiral phase transition at 
high baryon chemical potential from the continuous crossover at high temperature. Our 
attentions are paid to investigate the behaviors of the conserved charge fluctuations near 
the critical point, especial their singular behaviors, which shed light on the universal sym- 
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metry property of the QCD critical point. Furthermore, the singular behaviors of conserved 
charge fluctuations are also very helpful for searching for the QCD critical point in exper- 
iments [34|. In Fig. [T] we show the quadratic, cubic, and quartic fluctuations of baryon 



number as functions of \x& ({Iq = /is = 0) at several values of temperature calculated in the 
PNJL model. From the top panel of Fig. [TJ one can see that xi ^ as a P ea k structure when 
the chiral phase transition takes place, and this peak becomes sharper and narrower while 
moving toward the QCD critical point. We can also clearly notice that the xf diverges at 
the critical point, which is explicitly seen from the curve with the temperature of 160 MeV. 
According to the definition of cumulants of the conserved charge multiplicity distributions 
in Eq.(pQ), we have xi — dx% ' /d(fig/T). Since xi develops a cusp during the chiral phase 
transition, xi changes its sign there, which is clearly shown in the middle panel of Fig. [TJ 
This structure of xi was also found by Asakawa et al. 32j, who argued that the two sides 
of the QCD phase boundary can be distinguished by the sign of xi > therefore the third 
cumulants carry more information than the second ones. One can also find that xi diverges 
at the QCD critical point. Furthermore, we also calculate the fourth-order fluctuations of 
the baryon number, and the results are shown in the bottom panel of Fig. [TJ Based on the 
above analysis, it is expected that xi carries even more information than xi > since xi nas 
two positive maxima and one negative minimum as shown in Fig. [TJ Furthermore, we find 
that all amplitudes of xi i xii an d Xa grows rapidly when moving toward the QCD critical 
point and finally diverge there. 

In Fig. [2] we plot contours of the quadratic, cubic, and quartic fluctuations of the baryon 
number as functions of T and /zb calculated in the PNJL model. It is clearly seen that the 
chiral phase transition line in each of the three contour plots is obvious, and more important 
is that the region near around the QCD critical point also becomes manifest in the three 
plots, where the contour lines are dense. Therefore, our calculations demonstrate that by 
employing the second and higher order cumulants of baryon multiplicity distributions, it 
is possible to search for the QCD critical point. We should emphasize that in heavy ion 



collision experiments finite size and time effects should also be included [28|, |34( . Comparing 
higher order cumulants xi an d Xa with the quadratic one xi 1 we observe that the former are 
superior to the latter in the search for the QCD critical point. This is because only when the 
location is very near the chiral phase transition line, xi an d Xa are nonvanishing, while xi 
still has finite value when the location is far from the chiral phase transition line and in the 
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FIG. 2: (color online). Contour plots of quadratic (top), cubic (middle), and quartic (bottom) 
fluctuations of the baryon number as functions of temperature T and baryon chemical potential 
Hb (fJ>Q = fJ>s = 0) i n the PNJL model. 



chiral symmetry restored phase as shown in the top panel of Fig. El We should emphasize 
that in the top panel of Fig. HI there is a gray region in the lower right corner, which is 
because the quadratic baryon number fluctuation in this region have relative large value. 
However, this region does not correspond to any singular behavior, since \2 m this region 
is not divergent and it has a weak dependence on the temperature and baryon chemical 
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potential. This characteristic is quite different from that in the region near the QCD critical 
point, where xi changes rapidly with T and [Mb, and diverges at the QCD critical point, 
which is reflected by dense contour lines near the critical point in our plots. 
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FIG. 3: (color online). Ratios of cubic to quadratic (left panel) and those of quartic to quadratic 
(right panel) baryon number fluctuations as functions of baryon chemical potential [i B (pq = fj,g = 
0) at several values of temperature in the PNJL model. 

In Fig. [3] we show xi I X? an d xf I xf versus baryon chemical potential at several values 
of temperature calculated in the PNJL model. From Eq.(j3]) to Eq.(j5]) we find 

(SN B 3 ) 



xf 

Xf 

xf 

Xf 



2\ ' 



(SN B 2 ) 
(5N B 4 ) - 3(5N B 2 f 



(16) 
(17) 



(SN B 2 ) 

In fact, xi/xi anc ^ xf/xf are the skewness and kurtosis of the baryon multiplicity dis- 
tributions, respectively, which can be extracted from event-by-event fluctuations in heavy 
ion collision experiments. Furthermore, these ratios are intensive quantities, i.e., being in- 
dependent of the volume of the system. In Fig. [3] we show that xi/xi an d xf Vxf have 
the same structure as xf an d Xa > respectively, i.e., xf/xf change its sign during the chiral 
phase transition and there are two positive maxima and one negative minimum on xf/xf- 
Furthermore, it is observed that the amplitudes of xf/xf an d xf/xf grow rapidly when 
moving toward the QCD critical point and become divergent at the critical point. 

Fig. H] shows the quadratic, cubic, and quartic cumulants of the electric charge multi- 
plicity distributions as functions of fx B calculated in the PNJL model. It is seen that X2 
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FIG. 4: (color online). Quadratic (top), cubic (middle), and quartic (bottom) fluctuations of 
electric charge as functions of baryon chemical potential fis (a*q = Us = 0) with several values of 
temperature in the PNJL model. 

increases with the baryon chemical potential and has a cusp during the chiral phase transi- 
tion. Comparing X2 with xf > we observe that the peak in X2 grows less rapidly than that 
in X2 ■ A s 101 higher order fluctuations of the electric charge, it is found that xi changes its 
sign during the chiral phase transition and xl has two maxima and one minimum, which 
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are similar to xi an d Xai respectively. We should emphasize that once the location deviates 
from the chiral phase transition line, the higher order fluctuations of electric charge approach 
zero rapidly, while the second-order fluctuation X2 s ^ nas finite value even the location is 
far away from the chiral phase transition line. 




^ B (MeV) 

FIG. 5: (color online). Contour plots of quadratic (top), cubic (middle), and quartic (bottom) 
fluctuations of electric charge as functions of temperature T and baryon chemical potential \xb 
(hq = us = 0) in the PNJL model. 

In Fig. we show the contour plots of the quadratic, cubic, and quartic fluctuations of 
electric charge as functions of temperature and baryon chemical potential calculated in the 
PNJL model. From the three plots one can easily recognize the chiral phase transition line, 
which shows the same features as the contour plots of baryon number fluctuations given 
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in Fig. |2j However, we find that employing the quadratic fluctuations of electric charge to 
search for the QCD critical point is not easy, since the critical point in the top panel of 
Fig. is not obvious. While for the higher order fluctuations of electric charge, it is seen 
that the QCD critical point in the contour plots of \i an d xl is quite distinct. Therefore, 
higher order fluctuations of electric charge are more appropriate for being used to search for 
the QCD critical point in heavy ion collision experiments. 

Fig. [6] shows the quadratic, cubic, and quartic fluctuations of strangeness versus the 
baryon chemical potential at several values of temperature calculated in the PNJL model. 
It is noticed that the second-order and higher order fluctuations of strangeness are also 
enhanced when moving toward the QCD critical point, which are the same as the fluctu- 
ations of baryon number and electric charge. However, contributions to the singularity of 
the strangeness fluctuations from the QCD critical point are much less than those to the 
singularity of the baryon number or electric charge fluctuations. In the same way, it is found 
that higher order fluctuations of strangeness are superior to the second-order one in search 
for the QCD critical point. 



IV. NUMERICAL RESULTS OF CORRELATIONS OF CONSERVED CHARGES 

In this section, we are considering the correlations between or among conserved charges 
in the 2+1 flavor PNJL model. We will show that higher order correlations of conserved 
charges are sensitive to the critical behaviors related to the QCD critical point and therefore 
are very appropriate for being employed to search for the critical point. In Fig. [7] we show 
the dependence of the second-order correlations, i.e., Xii> X11 \ an d X11 011 the baryon 
chemical potential at several values of temperature in the PNJL model. It is found that 
all the second-order correlations have a non-monotonic behavior as functions of the baryon 
chemical potential. There is a peak structure on Xxi during the chiral phase transition; Xii 
is negative in the whole chemical potential region and has a minimum at the phase transition; 
and Xn also has a negative minimum but its value is positive in the chiral symmetry restored 
phase. Comparing Xii with X\\ an d X11 one can see that Xx\ does not vanish only when 
the thermodynamical system is very near the chiral phase transition, while the magnitude 
of Xii an d X11 increases with the baryon chemical potential in the chiral symmetric phase 
as shown in the middle and bottom panels of Fig. [7J This is because in the chiral symmetric 
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FIG. 6: (color online). Quadratic (top), cubic (middle), and quartic (bottom) fluctuations of 
strangeness as functions of baryon chemical potential \xb {^q = Us — 0) with several values of 
temperature in the PNJL model. 



phase, the system can approximately be described as noninteracting massless gases, i.e., the 
Stefan-Boltzmann limit. It can easily be shown that in the Stefan-Boltzmann limit Xii is 



vanishing while xfi an d Xii have finite values 



Fig. E] shows the contour plots of the second-order correlations Xii > Xii ' ■> an d Xii as 



as 
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FIG. 7: (color online). Second-order correlations between baryon number and electric charge (top), 
baryon number and strangeness (middle), electric charge and strangeness (bottom) as functions of 
the baryon chemical potential [ib (fJ-Q = Us = 0) with several values of temperature in the PNJL 
model. 

functions of temperature and baryon chemical potential calculated in the PNJL model. It is 
found that the chiral phase transition line in the three plots is distinct, but the QCD critical 



point in the plot of xfi is rnore apparent than that in xfi or Xii ■ 
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FIG. 8: (color online). Contour plots of the second-order correlations Xii (top), xfi (middle), 
and X\i (bottom) as functions of temperature T and baryon chemical potential \xb (mq = Us = 0) 
in the PNJL model. 

In Fig. we plot the third-order correlations X21 > X12 ' 1 X21 1 X12 ^ X21 > anc ^ X12 as 
functions of the baryon chemical potential at several values of temperature calculated in the 
PNJL model. It is seen that all these third-order correlations change their signs at the chiral 
phase transition, which are the same as the third-order fluctuations of conserved charges. 
More concretely, X21 an d X21 change their signs from negative to positive with the increase 
of the baryon chemical potential during the chiral phase transition, while other correlations 
in Fig. [9] changes their signs in the opposite direction. Furthermore, one can observe that 
the oscillating amplitudes of these third-order correlations all increase rapidly when moving 
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FIG. 9: (color online). Third-order correlations X21 (top-left), X12 (top-right), X21 ( m iddle-left), 
X\2 S (middle-right), X21 (bottom-left), and X12 (bottom-right) as functions of the baryon chemical 
potential \xb (hq = = 0) with several values of temperature in the PNJL model. 

toward the QCD critical point and diverge at the critical point. 

Fig. [TO] shows the contour plots of the third-order correlations X21 \ X12 \ X21 ' 1 X12 > 
X21 > an d X12 as functions of temperature and baryon chemical potential calculated in the 
PNJL model. Comparing Fig. [TU] and Fig. [HI one can easily notice that the QCD critical 
point in Fig. [10] is much more obvious than that in Fig. [8] which means that higher order 
correlations are more sensitive to the critical behavior of the QCD critical point than the 
quadratic correlations, and are better to be used for exploring the critical behavior of the 
QCD critical point in heavy ion collision experiments. As for the six contour plots in Fig. [JTJJ, 
it is easily seen that the QCD critical point in X21 > X21 ' 1 anc ^ X21 ^ s more distinct than that 
in Xi2 Q , Xii, and x? 2 S 
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FIG. 10: (color online). Contour plots of the third-order correlations X21 (top-left), X12 (top- 
right), xfi 5 (middle-left), X12 (middle-right), X21 (bottom-left), and X12 (bottom-right) as func- 
tions of temperature T and baryon chemical potential hb (mq = fJ-s = 0) in the PNJL model. 

In Fig. [TTIwe present the last third-order correlation Xiu_ S \ i- e -> the correlation among the 
baryon number, electric charge, and the strangeness, as a function of the baryon chemical 
potential with several values of temperature calculated in the PNJL model. We see that Xin S 
changes its sign from negative to positive during the chiral phase transition and diverges 
at the QCD critical point. Furthermore, it is noticed that only when the thermodynamical 
system is near the chiral phase transition, Xin has nonvanishing value. We also show the 
corresponding contour plot of the Xin as the function of T and in Fig. [T2J and it is seen 
that the QCD critical point in the contour plot of Xin S i s ver y obvious. Therefore, xfif is 
an ideal probe to search for the QCD critical point in heavy ion collision experiments. 
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FIG. 11: (color online). Third-order correlation among baryon number, electric charge, and 
strangeness, i.e., Xiu. S as a function of the baryon chemical potential [ib (fJ-Q = fJ-s = 0) with 
several values of temperature in the PNJL model. 
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FIG. 12: (color online). Contour plot of the third-order correlation Xui as ^ ne function of 
temperature T and baryon chemical potential fis (hq = f^s = 0) m the PNJL model. 
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V. SUMMARY AND DISCUSSIONS 



We have studied the fluctuations and correlations of conserved charges, i.e., the baryon 
number, the electric charge and the strangeness, in the 2+1 flavor Polyakov-Nambu-Jona- 
Lasinio model at finite temperature with nonzero baryon chemical potential. More attentions 
have been paid on the studies of the non-monotonic behavior of the fluctuations and cor- 
relations of conserved charges near the QCD critical point. The fluctuations are calculated 
up to the fourth-order and the correlations to the third-order. 

It has been shown that the second-order fluctuations and correlations have a peak or 
valley structure when the chiral phase transition takes place with the increase of the baryon 
chemical potential. As for the higher order fluctuations and correlations of conserved charges, 
it has been seen that the third-order fluctuations and correlations change their signs during 
the chiral phase transition and the fourth-order fluctuations have two maximum and one 
minimum. Furthermore, it has been noticed that the absolute values of the extrema of the 
fluctuations and correlations at the chiral phase transition increase rapidly when the thermo- 
dynamical system moves toward the QCD critical point (in heavy ion collision experiments, 
which means that the freeze-out point moves toward the QCD critical point), and finally 
the fluctuations and correlations of conserved charges diverge at the critical point. 

We have also explicitly demonstrated the critical behavior by depicting contour plots 
of the fluctuations and correlations of conserved charges. In these contour plots, one can 
clearly figure out the chiral phase transition line. Comparing with the second-order fluctu- 
ations and correlations, we have found that higher order cumulants, such as the third- and 
fourth-order fluctuations and the third-order correlations discussed in this paper, are more 
sensitive to the critical behavior of the QCD critical point. Therefore, we arrive at the con- 
clusion that the higher order fluctuations and correlations of conserved charges are superior 
to the second-order ones to be used to search for the critical point in heavy ion collision 
experiments. Particularly, we would like to address that among all the fluctuations and 
correlations discussed in this paper, the numerical calculations within the 2+1 flavor PNJL 
model indicate that X21 ' 1 X21 ' ■> X21 ' ■> anc ^ Xiu S are the most valuable probes for exploring 
the QCD critical point. 
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